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It was recently suggested that the error with respect to experimental data in nuclear mass cal- 
culations is due to the presence of chaotic motion. The theory was tested by analyzing the typical 
error size. A more sensitive quantity, the correlations of the mass error between neighboring nuclei, 
is studied here. The results provide further support to this physical interpretation. 
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The precision of nuclear mass spectrometry has dra- 
matically improved in recent years. Thanks to the Pen- 
ning trap technique, mass measurements with relative 
precision as high as 10~*, or even 10^^" for stable nuclei, 
may be achieved 0. This unprecedented accuracy has 
important consequences in different areas, like the deter- 
mination of fundamental constants, symmetry violations, 
metrology, stellar evolution and nucleosynthesis. 

An additional motivation for such precise measure- 
ments is to establish accurate and predictive mass for- 
mulae. Global nuclear mass calculations have been pur- 
sued over the years with increasing accuracy (see 
Ref.[3 for a review on recent experimental and theoret- 
ical developments). Despite the numerous parameters 
contained in the different models and the variety of the 
approaches adopted, two peculiar features emerge from 
these calculations. On the one hand, different models 
yield similar results for the known masses. A typical ac- 
curacy is 5 X lO^"' for a medium-heavy nucleus whose 
total (binding) energy is of the order of 1000 MeV. On 
the other hand, the predictions of different mass models 
strongly diverge when applied to unknown regions (they 
may differ by several MeV, i.e. relative variations of order 
5 X 10-3). 

These two features point towards the possibility of a 
basic underlying physical mechanism not appropriately 
incorporated into the present models. This mechanism 
should explain, in particular, the observed differences be- 
tween measured and calculated masses. In Ref.'Sl it was 
shown that the presence of chaotic layers in the nucleonic 
motion, whatever its physical origin may be, leads to a 
contribution to the nuclear mass whose typical size (Jch 
is given by 
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where A is the mass number. Equation J^l, obtained 
through a mean-field theory, follows from very general 
arguments, and is independent of any detailed informa- 
tion concerning the system. It provides an order-of- 
magnitude estimate of the chaotic contribution, and de- 
termines the onset of a new regime. Figure 1 shows, as 



FIG. 1: Root mean square of the difference between experi- 
mentally measured and calculated masses. Dots, squares and 
crosses from the calculations of Refs.0, and Q, respec- 
tively. The solid line shows the chaotic contribution, Eq.Q. 



a function of A, the comparison between this prediction 
and the typical size of the error of several nuclear mass 
formulae using the most recent experimental compilation 
Three different global mass models are compared. 
Two of them are based on mean field theory. The first one 
is a non-self-consistent macroscopic-microscopic model 
0, the second one is a self-consistent calculation based 
on Hartree-Fock-BCS j3| while the third one is a shell- 
model based calculation with parameterized monopole 
and multipole terms. The agreement with Eq.l^ for the 
two mean-field models is remarkable. The mass num- 
ber dependence of the error is also well described for the 
third model with, however, a factor of order two between 
their amplitudes, a tendency that one might expect for 
a model that includes substantial residual interaction ef- 
fects. 

When the difference between calculated and measured 
masses is plotted as a function of the neutron number (or 
of any other relevant parameter), an oscillatory curve, 
whose typical amplitude is shown in Fig.l, is observed. 
This curve doesn't look as a random white noise signal 
but shows structures. This qualitative remark is consis- 
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tent with theoretical expectations. Indeed, it is known 
from semiclassical mean-field theories |^ that the fluc- 
tuations of the mass produced by chaotic layers are dom- 
inated by the short classical orbits. The statistical prop- 
erties of these orbits show system dependent features. As 
a consequence, definite non-random fluctuations for the 
chaotic component of the mass are expected. Some of 
the predictions, like the asymmetry and non-Gaussian 
nature of the probability distribution of the oscillations 
0, have been tested recently j^. Thus, in the present 
context chaos should not be assimilated to a random un- 
predictable process. Though it may be a difficult task to 
explicitly compute its contribution for each nucleus \u\ 
because, generically, the detailed oscillatory structure of 
this contribution is very sensitive to details of the Hamil- 
tonian, it does not set an a priori bound to the accuracy 
of the theoretical mass calculations j^. This view dif- 
fers from that expressed in Refs. [l^ . where it is stated 
that the (partially) chaotic nuclear motion imposes limits 
in principle on the accuracy of the calculations of bind- 
mg energies (cf Ref.[i3). Moreover, and in contrast to 
global microscopic mass models, seemingly random devi- 
ations of typical size ~ 100 keV are observed in algebraic 
mass relations designed to locally cancel the interactions 
between nucleons 

Our aim here is to further develop the analysis of the 
mass errors in terms of chaotic motion by computing the 
autocorrelation function at different neutron numbers. 
Compared to the estimate of the typical error size, this 
is a finer analysis of the fluctuations, that tests more 
subtle dynamical information. As discussed below, by 
assuming the errors are due to the contribution of chaotic 
motion it is possible to obtain definite predictions which, 
in particular, fix a typical scale of the correlations. 

Fluctuations of the nuclear masses may be written, in 
a semiclassical mean-field expansion, as 
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The sum is over all the primitive periodic orbits p (and 
their repetitions r) of a classical underlying effective 
single-particle Hamiltonian. Each orbit is character- 
ized by its action Sp, stability amplitude Ap ^., period 
Tp = dSp/dE and Maslov index Vp ^.- x is a parameter 
on which the effective potential depends on. Though we 
let it for the moment unspecified, it will be chosen be- 
low to be the number of neutrons. The orbits entering 
this expression are all evaluated at the Fermi energy Ep . 
The latter is related to the mass number through the 
condition J^^^ p(£', x) dE — A^ where p is the average 
single-particle density of states. 

When the parameter x is varied, the correlation func- 
tion of the energy fluctuations U is defined as, 



where the brackets denote an average over a suitable pa- 
rameter window, which is large compared to the typical 
scales of oscillation of U and small on a classical scale. 
The correlation is evaluated using Eq.(|2Il. When the 
parameter x varies, and for a large number of particles 
where Sp 3> ?i, the main contributions to the variations 
of U come from the variations of the action. In a linear 
approximation, valid when the variations are large com- 
pared to h but small compared to S'p, the action varies 
as S{xq±x/2) ~ S{xq)±QpX /2, where Qp = dSp/dx\xo- 
Then, from ^ and l^l, dropping terms whose average is 
zero, we obtain 
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The double sum contains interfering terms between dif- 
ferent orbits. However, as shown in Ref. 8], it is not 
necessary to consider them. The main contributing or- 
bits are in fact the shortest ones (i.e, those having the 
shortest period). Their contribution is well approximated 
by taking into account only the diagonal terms p = p' , 
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In Eq.® we have used the semiclassical expression of 
the spectral form factor, Eq.® below, to express the 
autocorrelation in terms of an integral over time. The 
average over the cosine function is computed over the 
whole set of periodic orbits of period between r and r -t- 
dr. Denoting Pt{Q) the distribution of the values of Q 
of all these orbits, then 



C(x) = {U{xo - x/2) U{xo + x/2))x, , 



(3) 



(cos(Qa;/;i))^ = / cos {Q x / h) Pr{Q)dQ. (7) 

Equations lO and © are valid for both regular and 
chaotic motion. To proceed further, we must specify the 
type of dynamics we are considering. As in 0, we iden- 
tify the error between experimental (Ucxp) and calculated 
mass (Locale) as originating from the chaotic contribution 
([/) of the nuclear motion, Ucxp — C^caic + U. When the 
motion is chaotic, the number of periodic orbits having 
period between r and r + grows exponentially with r. 
There are numerical as well as theoretical evidences indi- 
cating that the distribution function Pt{Q) is in this case 
Gaussian, of average {Q)t = rjT and width {Q'^)t — olt 
jl4| . Nontrivial dynamical information is contained in a. 
When X represents variations of the neutron number, the 
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parameter rj takes into account the effect of the increase 
of the volume of the nucleus as a neutron is added, im- 
plying an increase of the average action (or length) of the 
orbits. Finally, an explicit expression for K£,{t) is also 
needed. The semiclassical expression for Kd{t) contains 
detailed information about the discrete spectrum of pe- 
riods Tp of the periodic orbits. 



(8) 
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Here, we will not take into account this detailed informa- 
tion, but instead use, as in 6], a continuous approxima- 
tion 
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where Tmin is the period of the shortest periodic orbit of 
the system. This schematic approximation incorporates 
two important features. On the one hand, it contains a 
system dependent information, namely the truncation for 
times smaller than Tmin- On the other hand, it displays 
the universality observed in chaotic systems with time 
reversal invariance for times Tmin ^ t ^ /ip, namely the 
linear growth characteristic of random matrix theory. 

Using the Gaussian form of Pt{Q) in Eq.jT)) as well as 
the approximation of Ko (t) , the correlation function 
Eq.(jni, normalized to one at the origin, takes the form 
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where the dimensionless parameters ( and a are given by 
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Following the lines developed in Ref.Q, a computa- 
tion of the average increase of the orbit's action when 
the volume V of the nucleus changes by 5V gives 77 = 
Ep5V/V. Computing 5V when a neutron is added, 
and using the expression Tmin/ot — (ttct^:)"^ |^, where 
cr^ = {{dUldxf)^'^, we obtain a = (\/2/7rA) E^/a^. A 
reasonable estimate for Ux, that has been tested numer- 
ically, is ~ 3/^1/3 MeV. Then a - / {2ttA^/^), which 
is of order 0.5 for A w 50 and 0.2 for A w 200. We have 
verified that, for these values of a, the error introduced 
by setting a = in Ea. llU|) is less than 7 x 10"^ for any 
C. Setting a to zero, the final result can be expressed as 



C(C) - ( 1 - ^ ) 
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where r(s, z) — t^~^e~*dt. Through the 

reparametrization all the system specific features 

have been incorporated in C. This leads to a "universal" 



function, Ea. H12() . Why it is so, as well as the validity of 
this result, is discussed below |l5l |. 

The parameter a contains detailed physical informa- 
tion related to the single-particle spectrum How- 
ever, it is difficult to extract the relevant information 
from experimental data. To circumvent this difficulty, 
the alternative expression Q 
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where OxqU — dU{xo)/dxo, is more convenient for 
our purpose because it only involves properties of U. 
Though the structure of Ea. (|13() is reminiscent of the 
reparametrization introduced in the context of random 
matrices and universal parametric correlations |l4l , 
there are however important differences: in the latter 
case single-particle energies are considered and universal- 
ities are exhibited, whereas here parametric correlations 
of thermodynamic properties of a Fermi gas are consid- 
ered, and important system specific features are shown 
to emerge (see below). 

From data, we have analyzed the following correlation 
function 
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where U{Z, N) is the difference between calculated and 
measured masses for a nucleus having Z protons and N 
neutrons, and dN is the difference in neutron number 
along an isotopic chain (a similar analysis for isotonic 
chains can also be performed). For a specific dN and Z 
every available f{N) = U{Z,N)U{Z,N + dN) is calcu- 
lated. To compute the mean value with respect to N one 
has to sum all the /(A^)'s and divide by the total number 
of nuclei in the isotope chain. This normalization ensures 
the non-negativity of the Fourier transform of the cor- 
relation function. The average (u'^'^ is obtained by 
setting dN = 0. In order to improve the statistics, that 
are severely limited by the length of the isotopic chains, 
a further average is computed over neighboring chains. 



C{dN) = {CzidN)) z. 



(15) 



Finally, to compare with Ea. H12|l . the correlation C{dN) 
should be expressed in terms of the parameter (. For each 
dN, (locally) the value of ^ is computed from Ea. lfT^ . 
where dxgU — U{Z, N +1) — U{Z, N), and averages over 
N and Z are as before. 

Figure 2 shows the normalized experimental correla- 
tion function of the nuclear mass errors, Ea. H15|l . plotted 
as a function of C, compared to the prediction (|12|l . The 
three different models studied give very similar results, in 
good agreement with theory. The overall quality of the 
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FIG. 2: Correlation of the mass errors between neighboring 
isotopes as a function of the dimensionless parameter Cir- 
cles, squares and crosses correspond to the same models as in 
Fig. 1. The solid line is the theoretical result, Eg. (1121 . 

prediction is comparable to Fig.l. The independence of 
the result with respect to the model used is a strong sup- 
port for our interpretation, though the influence of effects 
beyond mean-field theory still need to be clarified. 

Though Eq. H12|) is a continuous function of ^ (in par- 
ticular, it behaves as C(C) « 1 — close to the origin), 
we have chosen to plot it at discrete values of C, with a 
step similar to the experimental one. The half-width of 
the correlation function C{() is of order 1.5, which corre- 
sponds to neutron differences dN « 2, giving the parame- 
ter range over which the chaotic contributions of different 
isotopes are correlated. The existence of correlations ex- 
tending over a few nucleons seems to be consistent with 
the high accuracy with which Garvey-Kelson-type mass 
relations are fulfilled as mentioned before. 

The main approximation in Eq. (|12|l is the replacement 
of the diagonal form factor, Eq.ljSJ, by the function l^. 
By this, the system specific spectrum of periodic orbits is 
replaced by a function that keeps only one relevant scale, 
the period Tmin of the shortest one. It is this simplifica- 
tion that, through the reparametrization Hll|l . allows to 
obtain a "universal" correlation function C(C). Although 
this leads to a reasonable approximation, at least for the 
available neutron differences, the exact form of the cor- 
relation function of the chaotic contribution to the mass 
would be better described by the discrete sum (|SJ). This 
sum depends on the precise properties of the periodic or- 
bits, and is therefore system and model dependent. Be- 
cause of its discrete nature, generically oscillations of the 
correlation as a function of x are expected. This contrasts 
with the uniform "universal" non-oscillatory decay given 
by Eq. H12() , which is clearly an artifact of the approxima- 
tion @. Some tendency towards negative values, and 
therefore of oscillatory behavior, seems to be present in 
the autocorrelation of the errors shown in Fig. 2 at large 
values of However, to be conclusive, larger values of C, 
which are not experimentally available, are needed (this 



problem might be overcome by considering different nu- 
clear chains 

To conclude, let us recall that shell effects described by 
periodic orbit theory related to regular motion are very 
familiar in nuclear physics, particularly since the work of 
Strutinsky (see Ref.[l3 for a recent discussion). What 
about the contribution of unstable chaotic orbits? It has 
been recently suggested that the size of the differences 
between measured and computed binding energies can 
be attributed to the presence of nuclear chaotic motion. 
The work presented here goes one step further in this 
direction by studying the autocorrelation function of the 
chaotic contribution to nuclear masses. The result, that 
depends on some well identified physical parameters, is in 
good agreement with the autocorrelation computed from 
several existing models along isotopic chains. This gives 
further support to the view that chaotic dynamics effects 
are present in the ground state of atomic nuclei. 
* Unite de recherche associee au CNRS. 
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